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Appendix  I 


Historical  Remarks  on  HOPF  Bifurcation  in  Aircraft  Dynamics 


Modern  qualitative  theory  of  differential  equations  originated 
with  the  work  of  Poincare,  who  essentially  had  arrived  at  the 
phenomenon  which  is  presently  called  Hopf  bifurcation.  Poincare's 
work  was  later  developed  by  Androuov  and  Pontryagin  [l]-[3]]  and  Hopf 
[15]. 


More  recently,  experience  with  the  aircraft  at  high  angles  of 
attack  shows  loss  of  stability  for  some  critical  values  of  parameter 
(e.g.,  angle  of  attack,  or  the  velocity).  This  phenomenon  was 
interpreted  by  R.  Mehra  et  al.  [25]  as  the  Hopf  bifurcation;  this 
interpretation  was  supported  by  numerical  study  for  the  aircraft  H 
model . 


The  work  of  Mehra  et  al.  led  naturally  to  the  question  as  to  the 
effect  of  control  decoupling  feedback  on  the  (undesirable)  Hopf 
bifurcation;  in  particular,  does  such  a  feedback  eliminate  the 
bifurcation?  Among  other  things,  this  question  will  be  answered  in 
this  report. 


The  problem  of  control  decoupling  via  feedback  together  with  the 
closely  related  problems  of  disturbance  decoupling  and  Invariance  has 


arisen  In  many  engineering  applications,  and  particularly  In 
connection  with  the  aircraft  control  problem.  Since  Rozenoer's 
Initial  work  [25],  the  subject  of  control  decoupling  via  feedback  has 
been  extensively  developed,  and  a  reasonably  large  body  of  literature 
currently  exists.  Some  papers  which  have  been  important  milestones 
are  Wonham  and  Morse  [33],  Tokamaru  and  Iwal  [32],  Majumdar  and 
Chaudhury  [17],  Isidori  et  al.  [16],  Hischorn  [14],  and  Byrnes  and 
Krener  [7].  In  addition  to  this  theoretical  work,  several 
applications  to  prohleras  of  aircraft  control  have  been  studied 
Including  Singh  and  Schy  [27]  and  work  by  G.  Meyer  [20]  on  the  design 
of  an  autopilot  system  for  the  Bell  UH-IH  helicopter. 

The  basic  Idea  behind  the  theory  of  control  decoupling  Is  quite 
simple:  Suppose  there  is  given  a  nonlinear  control  system  of  the  form 
x  -  f(x)  +  u^Cx)  +  u  2®  2<  x ) 
yi  "  H1(X)  (i“1.2) 

We  wish  to  consider  modifications  of  the  system  dynamics  using  feedback 
controls  u»  a(x)  +  8(x)v 


f(x)  +  viRl(x)  +  V282^X) 


x 


(2.1) 


7  "  h  (*)  (i-1,2)  (2.2) 

where 


2 

f  -  f  +  E  g. 

j=X  3  3 


2 

■  ;■ 

The  decoupling  problem  Is  to  find  a  and  8  such  that  controls  and 
only  y^.  (That  Is,  we  want  to  have  no  influence  on  the  output  y2 
and  vice-versa). 


Techniques  for  finding  a  and  8  in  the  case  in  which  f  and  the 
h^'s  are  linear  and  the  g^'s  are  independent  of  x  are  well  known  and 
may  be  found  in  Uonham  [33].  For  nonlinear  systems,  considerably  less 
is  known  and  although  the  beginnings  of  the  theory  date  back  to  1962, 
many  aspects  remain  to  be  understood. 


In  1962,  L.  Rozonoer  [25]  obtained  conditions  necessary  for 
Invariance  (i.e.  independence  of  the  output  upon  the  input)  in 
nonlinear  systems,  by  using  variational  method,  similar  to  the 
Pontryagln's  optimality  principle. 


3 


H.  Tokamaru  and  Z.  Iwai  [32],  and  A.  Najumdar  and  A.  Chaudry 
[17],  independently  have  applied  the  variational  method  used  by 
Rozonoer  to  obtain  necessary  conditions  for  non-interacting  control, 
for  linear  time-variant  systems  (1968)  and  for  nonlinear  systems 
(1971). 

A  concrete  design  for  control  decoupling  feedback  for  the 
automatic  piloting  of  a  mine  hunter  boat  was  given  by  E.  Daclln  [8]. 

All  theoretical  work  mentioned  above  neglected  the  existence  of 
constraints  on  controls  in  virtually  every  real-life  situation  (e.g., 
a  limit  on  rudder  deflections,  acceleration  rates,  etc.  The  effects 
of  such  constraints  have  been  studied  by  D.  Hanson  and  F.  Stengel 
[12],  particularly  for  systems  with  two  degrees  of  freedom.  Further 
development  of  their  analysis  seems  to  be  of  considerable  Interest. 

A  different  approach  to  the  problem  was  used  by  R.  Su,  G.  Meyer 
and  L.  Hunt  [31],  who  used  nonlinear  transformations  to  reduce  the 
nonlinear  problem  to  the  linear  one,  which  can  be  treated  by  standard 
methods,  as  in  Wonham  [33 j.  The  main  drawback  here  is  that  it  is  not 
generally  possible  to  carry  out  this  type  of  linearization;  there  are 
both  algebraic  and  topological  obstructions  which  are  essentially  the 
same  ones  that  are  encountered  if  one  were  to  try  to  do  nonlinear 


decoupling  directly 


An  important  application  of  control  decoupling  theory  to  the 
problem  of  aircraft  dynamics  was  given  by  S.  Singh  and  A.  Schy,  [27]. 
These  authors  derive  a  feedback  for  a  simplified  aircraft  model  in 
which  certain  aerodynamic  forces  had  been  Ignored.  One  of  the  alms  of 
this  report  is  to  produce  a  control  decoupling  feedback  for  a  full 
blown  model  -  the  so-called  aircraft  H  model.  In  Section  3.1  we  will 
give  what  we  believe  to  be  the  simplest  derivation  of  a  decoupling 
feedback  control  law  together  with  a  simple  geometrical  explanation  of 
the  method. 
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This  section  consists  of  three  parts.  In  the  first,  we 
Investigate  the  underlying  geometry  of  control  decoupling  feedback, 
with  the  aim  of  gaining  intuition  about  the  problem  of  Interest  -  the 
control  of  certain  model  of  aircraft  dynamics,  the  so  called  "aircraft 
H"  model. 

Second,  we  give  a  brief  outline  of  Poincare's  theory  of  normal 
forms,  which  is  used  to  study  the  Hopf  bifurcation.  We  follow  the 
work  of  Poincare,  Hopf,  Andronov  et  al.  and  Arnold. 

Third,  we  apply  the  above  two  ideas  to  study  a  model  with  three 
degrees  of  freedom,  to  illustrate  the  effect  of  control  decoupling 
upon  the  Hopf  bifurcation. 

Using  the  theory  of  normal  forms  our  goal  is  to  reduce  our  system 
to  a  simpler  form;  the  subsequent  application  of  control  decoupling 
feedback  will  illustrate  the  stabilizing  (or  destabilizing)  effect  of 
control  decoupling  upon  the  Hopf  bifurcation.  We  shall  also  derive  a 
criterion  for  determination  of  stability  of  the  decoupled  systems. 


3.1  Derivation  and  Geometry  of  a  Control  Decoupling  Feedback 

Here  we  provide  a  straight  forward  derivation  of  a  control 
decoupling  feedback  for  the  system 
x  -  F(x)  +  G(x)u 

y  -  Hx  (3.1.1) 

As  we  have  remarked  In  Section  2,  the  objective  Is  to  modify  the 
dynamics  using  state  feedback  in  such  a  way  that  the  1-th  control 
input  Influences  only  the  i-th  output. 

Differentiating  the  output  and  using  the  feedback  u»*  a+gv  (with  a(x), 
g(x)  to  be  found)  we  obtain 

y  -  Hx  -  HF  +  HGu  -  (HF+HGa)  +  HGgv 

Our  aim  Is  to  find  a  and  g  such  that  v^  controls  y^  and  only  y^. 
Such  a  choice  Is  determined  by  setting  y*Av  (A  Is  an  arbitrary 
constant  diagonal  matrix),  l.e.  HF+HGa=0,  HGg*=A. 


We  obtain: 
a  =  -  (HG)-1H  F 

(3.1.2) 

g  -  (hg)_1a 

the  desired  feedback  is:  u*-(HG)-^HF+(HG)"^av. 


3.2  Poincare's  Theory  of  Normal  Forms 


One  of  the  most  powerful  tools  in  studying  dynamical  system  is 
the  method  of  normal  forms,  which  Poincare  introduced  in  his 
dissertation.  This  method  will  be  briefly  described,  and  we  refer  the 
interested  reader  to  Arnold  (1983)  for  further  details.  The 
development  of  normal  forms  for  nonlinear  control  systems  can  provide 
considerable  insight  into  problems  of  synthesizing  stabilizing 
feedback  control  laws.  In  general,  by  a  "normal  form"  for  a  control 
system  we  mean  a  canonical  form  with  respect  to  local  dif feomorphic 
changes  of  state  coordinates  near  an  equilibrium  for  the  force-free 
"drift"  term.  In  terms  of  our  general  model,  this  drift  term  is  given 
by  the  vector  field  F(x).  The  main  idea  behind  Poincare's  theory  of 


normal  forms  is  to  reduce  a  system 

x  -  Ax  +  higher  order  terms  (x  eR  ) 


(3.2.1) 


to  as  simple  a  form  as  possible  by  means  of  transformations 
x*T(z)»  z  +  h.o.t.  Ideally,  one  would  like  to  linearize  (3.2.1),  but 
the  existence  of  such  a  linearizing  transformation  depends  upon 
certain  non-resonance  conditions.  That  is,  there  may  be  Diophantine 
relations  between  the  eigenvalues  of  A  which  prevent  one  from 
eliminating  some  higher  order  terms. 


More  precisely,  a  collection  of  eigenvalues  u=(u,  , . . • ,u  )  of  the 

1  n 

nxn  matrix  A  is  called  resonant ,  if  there  is  a  relationship  of  the 

form  u  “(m,u) ,  where  ra=(m,  , . . . ,m  )  with  m  >  0 ,  £  m,  >  2 .  Here  ( . , .) 
s  In  k—  K— 

denotes  the  usual  dot  product.  Such  a  relationship  is  called  a 


resonance.  For  example,  y^  +  y2*0  *s  a  resonance»  or  rather,  it 

Implies  the  resonance  y^=2y^  +  y2 »  and  more  generally,  yj*(k+l)y^+ 

ky2»  y2  -  ky^  +  (k+l)y2,  k-1,2 .  Let  x^...,  x^  be  the  coordinates 

In  the  eigenbasls  e  e  of  the  matrix  A.  Equivalently,  we  can 

1  n 

perform  a  preliminary  linear  transformation  diagonalizing  A;  then  eg 
will  have  coordinates 

e  -  col  (0,...,1  ,...,0),  and  A  =  diag  (y,,...,y  ). 
s  s  in 

Therefore,  there  is  no  loss  of  generality  in  assuming  that  A  Is 
diagonal  (we  omit  the  case  of  Jordan  blocks,  which  can  be  treated 
without  difficulty,  too). 

ml  mn 

Consider  a  vector  monomial  xme  ,  where  xm  =x,  . . .x  is  a 

s  in 

scalar  monomial,  and  e  is  an  eigenvector  of  A.  It  is  called  resonant 

s 

if  y  »(m,  y)  ,  m  >  2  . 

s  — 

Now,  the  main  theorem  (of  Poincare  and  Dulac)  asserts  that  the 
system 

x  «  Ax  +  . . . , 

can  he  reduced  by  a  (formal)  change  of  variables  x=y  +  ...  to  the  form 
y  *»  Ay  +  w(  y) 

where  all  the  vector  monomials  in  the  series  w  are  resonant.  We  give 
a  short  proof  of  this  theorem  in  Appendix  1. 
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We  Illustrate  Poincare-Dulac's  theorem  by  applying  it  to  the 
proof  of  the  existence  of  a  Hopf  bifurcation. 

Consider  the  systems  (cf.  Section  3.3) 
x  ■  Xx  +  wy  +  ... 

(3.2.2) 

y  ■  -wx  +  X  y  +  . . . 

where  X  is  a  parameter  near  0.  When  X»0,  we  have  a  resonance,  p  + 
“  0  2“  +  and  therefore  for  X*0  we  cannot  kill  the  resonant 

terms  of  third,  fifth,  etc.  order.  Therefore,  we  do  not  try  to  kill 
them  at  all,  even  for  X=0  -  an  attempt  to  do  so  would  lead  to  the 
discontinuous  dependence  of  the  vector  field  on  the  parameter  X. 
Eliminating  the  second  and  fourth  order  terms  by  a  change  of 
variables,  we  can  reduce  (3.2.2)  to  the  form 

*1  “  Axi  +  wyi  +  ^axi  +  6x2^xi  +  >  +  °5 

(3.2.3) 

yi  “  -wxi  +  Ayi  +  ^'Bxi  +  aX2^xx  +  )  +  o5 

This  reduction  also  justifies  our  seemingly  specific  choice  of 
the  example  in  Section  3.3. 

To  eludicate  the  geometrical  meaning  of  this  analysis,  we  refer 


to  Figure  3.2-1 


2  2 

Introducing  the  distance  r  *  (x^  +  )  from  the  origin,  we 

obtain  the  the  following  ODE  for  r: 
r  -  x  r  +  ar2+  0^ 

3 

where  0^  refers  to  the  terms  of  order  r  and  higher. 

Ignoring  0^,  we  obtain  the  equation 
r  »  (  X  +a  r)r  (r>0) 

for  which  r  q  *  0  Is  always  a  stationary  point  (which  corresponds  to 
the  equilibrium  solution  of  (3.2.2));  another  stationary  point 
r^  -  -Va»  (if  <0)  corresponds  to  a  periodic  solution  of  the 
truncated  system  (3.2.3)).  Actually,  the  presence  of  0,.  in  (3.2.3) 
will  not  destroy  this  solution  but  will  only  perturb  It  slightly,  as 
can  be  easily  shown.  This  periodic  solution  is  close  to  a  circle  of 
radius  £  V  a  ;  It  bifurcates  off  the  origin  as  \  crosses  0.  This  Is 
exactly  the  Hopf  bifurcation. 

During  almost  a  century  of  its  existence,  Poincare's  method  has 
been  used  by  many  researchers,  Including  Birkhoff,  Hopf,  Moser  and 
others,  In  conjunction  with  their  work  on  problems  arising  in  such 
applications  as  electrical  engineering  and  classical  and  celestial 
mechanics,  to  name  but  a  few. 


It  is  known  that  the  aircraft  H  model  of  section  4  enters  an 
oscillatory  regime  as  angle  of  attack  becomes  large.  In  Mehra,  Kessel 
and  Carroll  (1978)  it  was  demonstrated  that  this  instability  is  due  to 
Hopf  bifurcation.  Our  eventual  goal  is  to  determine  the  possibility 
of  stabilizing  this  system  via  decoupling  feedback.  (We  construct 
such  a  feedback  in  section  4,  by  applying  the  method  given  in  section 
3 . 1  above) . 

Before  attempting  this  derivation,  we  would  like  to  understand 
the  geometry  of  the  problem;  to  that  end,  instead  of  the  aircraft 
model,  as  a  first  step  we  consider  a  simpler  three-dimensional  system 
with  two  controls,  which  retains  the  essential  instability  of  the 
aircraft  model,  but  is  easier  to  visualize  geometrically. 

The  role  of  equilibrium  angle  of  attack  (or  equivalently  elevator 
angle  6e)  will  be  played  by  a  certain  bifurcation  parameter  X  in  our 
example.  In  analogy  with  the  angle  of  attack  in  the  aircraft  model, 
as  X  increases,  our  system  will  pass  from  stable  steady  state  to  an 
oscillatory  regime. 

We  will  then  introduce  a  control  decoupling  feedback  into  our 
system  and  will  determine  its  effect  upon  stability.  The  importance 
of  considering  this  example  is  that  is  suggests  the  direction  of 
analysis  for  the  actual  aircraft  model. 


Thus,  we  consider  a  family  of  systems 


(3.3.1) 

(3.3.2) 


x  -  F(x,X,a)  +  G(x)u, 
y  -  Hx 

2  2  1 

with  x  eR  ,  ueR  and  H  ■  0  1  being  the  matrix  of  projection  of  the 

0  0 

vector  (x^^.x^)  on  the  (x^  ,x2 ) -plane. 

The  parameter  X  will  be  roughly  analogous  to  the  equilibrium 
angle-of-attack,  while  a  will  characterize  the  geometry  of  the  system 
and  will  be  referred  to  later. 


One  of  our  major  techniques  for  analyzing  stabilization  schemes 
for  C3.3.I)— (3.3.2)  is  to  put  this  control  system  into  a  normal  form, 
by  a  change  of  coordinates  in  the  state  variables  which  is  locally 
defined  in  a  neighborhood  of  an  equilibrium  point.  Note  that  such  a 
change  of  coordinates  does  not  change  the  time  trajectories  (y(t), 
u(t))  of  (3.3.1)— (3.3.2) ,  but  rather  enables  us  to  define  a  feedback 
strategy  for  (3.3.1)  in  the  context  of  the  simpler,  normal  form. 


We  first  describe  the  free  motion  of  (3.3.1)  in  a  geometric 


fashion.  To  be  specific,  consider  the  third  order  control  system 

x,  =  Xx„  +  wcx  +  wsx  -  x  (x2+  (cx  +  sx  )2) 

1  1  2  311  23 

+  g  u  +  g  u 
11  1  12  2 

X„  *  -  U)CX,  4-  (  XC^-  B2)  X  +  (x+l)csx  -  c(cx  +  SX  )(x2+(cx  +SX  )2) 

2  1  2  323123 


X  ■  -  U)SX^  +  (X+Dcsx^  +  ( Xs  -  c  )  x^  -  s(cx2  +  BX^Kx^+CcX^SXj)  ) 

+  83^  +  832u2 

where  G(x)«(g^(x))  is  assumed  to  be  of  full  rank  in  a  neighborhood  of 
0,  gj  is  a  fixed  real  quantity,  c«cosa,  s*sina,  end  x  and  a  are 
parameters.  Consider  the  free  motion  of  this  system  (with  controls 

ul*u2*°) 

•  2  2 

XI  *  *X1  +  u)cx2  +  t*,8X3”  x^(x3+  (cxji-  SX3)  ) 

2  2  2  2 

X 2  “  **  (1)0X3  +  (XC  -  8  )  X2  +  (X+DC8X3  -  c(cx2  +  SX3)(Xj+( CX2+SX3)  ) 

2  2  2  2 

X3  ■  -  0)8X3  +  (X+l)csx2  +  (X8  -  c  )  X3  -  s(cx2  +  8X3)(X3+(CX2+8X3)  ) 

(3.3.4) 

Standard  results  (see,  for  example  Marsden  and  McCracken,  [  ])  show 

that  as  X  (a  real  parameter)  increases  from  negative  values  to 
positive  ones,  a  limit  cycle  for  (3.3.4)  appears.  (See  figure  3.3-1) 
Indeed,  if  we  fix  the  parameter  a,  this  system  represents  a  normal 
form  for  third  order  families  of  systems  admitting  a  Hopf  bifurcation 
The  limit  cycle  is  confined  to  an  invariant  two  dimensional  manifold 
Wc  (which  for  X“0  is  called  the  center  manifold).  The  other  invariant 
manifold,  Ws,  is  stable  for  all  X,  snd  it  intersects  Wc  transverally. 

Note  that  there  is  a  plane  through  the  origin  perpendicular  to  Ws,  and 
the  motion  of  the  system  restricted  to  this  plane  is  given  by 
/x\  /  xx  +  ujy  -  x(x2+y2)  \ 


-u ix.  +  Xy  -  y(x2+y 2)/,  where  x,y  are  coordinates  in  that  plane 


We  point  out  that,  although  this  picture  may  seem  to  be  very  specific. 


It  accurately  represents  what  actually  occurs  in  the  aircraft  model, 
except  that  in  the  latter  dimWS«5  instead  of  1  and  dim  (WS)J'«3  Instead 
of  2. 


According  to  section  3.1,  the  feedback 


u  -  -  (HG)-1H  F(x,  A  ,ot)  +  (H  G)-1  A  v 


decouples  the  closed  loop  system.  The  resulting  system  is 


*i  '  V 


*2  -  X2v 


*3  *  F3  '  61  F  "  B2  F2  +  61X1  'l  +  B2B2  v2 * 


/  i  o\ 


where  G(H  G) 


-1 


0  1 


\h  62  / 


F2  I  and  B2  are  functions  of  x 
F„ 


More  explicitly,  we  have 


xi  vi 


X2  v2 

Ax1  +  Bx,  +  Cx,  +  N(x,,  X,,  x,)  +  B.X  vL  + 


2  2  2' 


where 


A  ■  -(jS  -  \  +  B2ojc 

B  -  (xH)cs  -  p  (jc  -82<Xc2-  s2) 

C  -  (As2-  c2)  -  B^u)S  -  p  (  A+l  )cs , 


and  N  includes  cubic  terms 


P r=  T'  TT  ».  7WT* \.V  ^  '.V  ■A^AT*.* r *•■'  «'■  ^  ■  . 


Now,  to  determine  stability  of  this  system,  we  keep  x^-x^^O  (by 
setting  v^v^O  and  making  initial  values  of  x^,x2  vanish);  the 
resulting  system  for  x^  is  x^"C(x,a)x^+N(0,0,x^);  it  is  stable 
precisely  when  C(x*a)<0« 

Thus,  we  obtain  a  stability  criterion;  If  C(0,q)<0,  then  Hopf 
bifurcation  is  “stabilized"  by  the  control  feedback  (3),  l»e.  as  long 
as  x^iX2  remain  small,  so  does  x^. 

Explicitly,  this  criterion  is: 

2 

-c  -  B^ws  -  32u>c8  <  0 

The  importance  of  this  criterion  is  that  it  points  to  the  fact  that 

either  situation  can  occur:  the  feedback  may  either  stabilize  or 

destablillze  the  bifurcation.  For  instance  if  we  have  Bj*62“W"1,  the 

criterion  becomes 
o 

cos  a  +  sin  a  +  cosa  sina  >  0; 

for  some  angles  a(e.g.  o^O)  this  inequality  holds,  and  for  others 
(e.g.  a«-  ^  )  it  doesn't.  In  each  specific  case  one  would  have  to 
make  calculations  to  determine  the  sign  of  C. 


When  this  approach  is  generalized  to  the  aircraft  model,  we 
obtain  a  similar  criterion;  instead  of  the  sign  of  the  coefficient  C, 
we  have  to  check  stability  of  a  certain  4x4  matrix  (corresponding  to 
the  uncontrolled  variables).  The  details  of  this  are  sketched  in 
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section  4  where  the  state  dynamics  for  x  ,  x  and  x  is  decoupled,  but 

14  5 

the  state  dynamics  for  x^,  x 
Flight  experience  with  advanced  fighters  such  as  F-14  shows  that 
certain  control  designs  e.g.  aileron-rudder  interconnect,  can  produce 
problems  in  other  flight  regimes.  We  believe  further  research  along 
the  lines  we  have  indicated  will  lead  to  analytical  methods  for 
solving  such  problems  by  examining  the  global  stability  properties  of 
given  feedback  laws. 


y  Xgand  x ^  is  further  complicated. 


A.  Application  to  the  Aircraft  H  Model 


In  this  section  we  apply  the  technique  described  In  section  3-1 
above  to  the  problem  of  control  decoupling  in  the  full  aircraft  H 
model  and  give  the  method  to  determine  stability  of  the  resulting 
decoupled  system. 


The  seven  degree  of  freedom  aircraft  model  (Hecker  and  Oprisiu, 

1976)  has  the  form: 

p  **  £.6  +  £  Aa6a  +  £  q  +  £  r  +  £  BAa  +  £  rAa  +  £  p  -  i.qr  +  £  6a  +  t,  Sr 

6  a  6a  q  r  Ba  ra  p  1  6a  6r 

q**mAa  +  mq+i„pr  +  m_  6e  -  m*pB 

a  q  2  6e  a 

r  *  n_B  +  n  .  Ao6a  +  nr  +  np  +  n  pAo  -  i_pq  +  n.  6a  +  n_  6r 
B  a  6a  r  p  par  6a  6r 

a  =  q  -  pB  +  z  Aa  +  z,  6e  +  (q/V)  (cos  0  cos  <J>  -  cos  0  )  (4.1) 

a  oe  o 

B  -  yaB  +  p(8in(a  +  Aa)  -  r  cos  a  +  y.  6a  +  y*  6r  +  (g/V)  cos  0  sin  4> 

p  o  03.  or 

41  “  P  +  q  tan  0  sin  <j>  +  r  tan  0  cos  <f> 

• 

0  *  q  cos  -  r  sin  41 

where  m  •»  m  +tn*z  ,  m  -m  +m*,  ra.  =m.  +ra»z.  ,  and  where  the  symbols  are 

a  a  a  a  q  q  a  6e  6e  a  6e 

a 8  follows: 


x,y,z  principal  axes  of  aircraft 

2 

I  ,  I  ,  I  moments  of  inertia  about  principal  axes  (kg/m  ) 


(i  -I  )/l 
z  y  x 


(I  -I  )/I 
z  x  y 


(I  -I  )/I 


•\ 


non-dimensional  inertia  coefficients 


p,q,r  Angular  velocity  components  along  roll,  pitch  and  yaw 


axes  (principal  axes)  of  aircraft  (rad/sec) 

V  velocity  of  aircraft's  center  of  mass  (m/sec) 

2 

g  gravitational  acceleration  (m/sec  ) 
a  angle  of  attack  (rad) 

3  angle  of  sideslip  (rad) 

0  pitch  angle  (rad) 

^  angle  of  bank  (rad) 


6a,  6 r. 


a  ,6  »a  »P  i<I  •  r 


fie  aileron,  rudder,  and  elevator  deflection  angles  (rad) 

o 

rolling  moment  per  I^(  l/sec  ) 

2 

m  pitching  moment  per  I  (l/sec  ) 

2 

n  yawing  moment  per  1^( l/sec  ) 

y  side  force  over  (aircraft  mass  times  speed)  (l/sec) 
z  aerodynamic  force  along  principal  axis  z  over  mass  times 
speed  (l/sec) 

,fia,fir,fie  as  subscripts  denote  partial  derivatives  with  respect  to 
the  respective  quantity,  i.e. 


2 

=  =  _  =  __9m _  n  _  3  n 

8  9  6  *  "6a  3  6a  ’a  ^  da  *  6aa  3a6a 

9  ■  j_— 

<t  bank  angle  command 

n 


We  change  the  notation  to  a  more  uniform  one  by  setting 
(p,q,r,a,  6,<j>,0)-  (xltx2,x3,x4,x5,x6,x7),  £p"£i^g  *l2'  etc*  (the 

index  of  each  coefficient  corresponds  to  the  index  of  the  variable  (or 
variables)  which  it  multiplies). 


Our  system  can  be  rewritten  as 

x  ■  F(x)  +  G(x)u  ,  y  **  Hx  (4.2) 


the  explicit  form  of  F^x)  is  evident  from  (4.1),  and  ai  are 
constants,  except  for  a^  and  a^: 


al 

■ 

i  c  Ax,  +  5-i_, 

a  6a  4  6z 

a2 

-  1. 

6r 

a3 

- 

6e 

a4 

■  n  „  Ax,  +  n 

a 6a  4  6a 

a5 

B 

n6r 

a6 

Z6e 

a7 

B 

y6a 

a8 

"  y6r 

v  _■  -•  .* 


V  >  ■ 


v.  -r. 


JljS 

Ciii 


Me 


W“,  .*, 

W' 


'.•-'S’ 

r-», 


=Wi 

K.'-v- 


W: 


t... 


m 


'- 

L»  *  -  ’  < 
F  ,  "  •  " 


Ijl  "  -  ^ 

kv->. 


.-  V 

A, 

►  V- 


*v\V 

r*  •  •  m  * 

k  -  -  •  - 

h#.'" 


;v\-r- 


v$i 


K-: 

K 


t’:::  :: 


The  procedure  described  in  section  3-1  leads  to  the  following 
results  (we  omit  the  tedious  calculations). 


The  functions  a(x)  and  $(x)  in  the  control  decoupling  feedback 
u  ”  a(x)+g(x)v  are  as  follows: 


\  0  h2 


B(x)  -  (HG)  -1A  -  A-1  0  h3  0 


h.  0  h, 
4  5/ 


»!*>!  0  x3h2  ^ 


0  xnh.  0 

J.  J 


Xlh4  0  X3h5/ 


(4  .3) ■ 


where  A  *  aia6ag'*a2a6a7  and  h^a^g,  h^-a^,  h^-a^+a^g,  h^-a^. 


h  **a  a 
5  16 


Vl+  V2 


a(x)  -  -( HG) -1H  F  -  A-1  F,h 


4  3 


(4.4) 


'h+Fh  / 
1  4  5  5  / 


The  feedback  u=a(x)+6(x)v  results  in  the  following  system: 


_d 

dt 


~X1 

0 

X1V1 

x2 

FrAiF4 

A1X2V2 

X  3 

F3_A2FrA3F5 

A2X1V1+A3X3V3 

X4 

= 

0 

+ 

X2V2 

X  5 

0 

V3 

X  6 

F  6 

0 

x7 

F7 

0 
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v.v  . . 


■  .  . 


■»v- 


>  ‘  .v  \  , 

•  *  «■  v'  •  *  v'  «y.  * 


A  -  v*.'  V  \’  «.*  d 

*  *  *  •  -  «L*  *  ’  fc  ‘V1'.  **« 


with  ajj3  standing  for  a^a^Sg,  etc»  Local  stability  of  the 


uncontrolled  part  of  the  system  depends  on  the  eigenvalues  of  the 


matrix  of  the  linearization  of  the  dynamical  system  restricted  to  the 


4-dimensional  subspace  given  by  x  =x  *0,  x  -x;  -2.5°  (the  leveled 
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steady  flight  angle  of  attack).  Explicitly,  this  dynamical  system  has 


the  following  form: 


Fr  aif4 


+  *1  X2V2 


F  -A  F  -  A  F  +  AX  v  +  A,  v 


3  2  1  3  5 


2  1  1  3  3  3 


Here  the  components  F  of  the  vectorfield  are  as  follows: 


F1  '  *5X5+  e2V  *3  V  *45nx5  ‘Vs'V  *lV  11X2X3 


F3  '  "sV  n3*3+  "iV  "l4VV  13X1X2 


F  ■  X  -  X  X  +  z,  &x,+  &  (cos  x,  cos  x,  -  cos  x°) 
421544V  7  6  7 


Fc  -  ycxc  +  x,  sin  x^,-  x^cos  x°  +  —  cos  x7  sin  x 


5  '5  5  1 


Fg  -  x^  +  x2  tan  x^  sin  x^  +  x^  tan  x^  cos  x^ 


F^  ■  x2  cos  Xg  -  x^  sin  Xg 


The  resulting  nonlinear  system  becomes,  after  setting  x^**Xg*0, 


x^-  x^( so  that  Ax^-0): 


x2  -  U2-  Aj)  x2  +  ^x3  -  i^x.,  ~y~  cos  x,.  cos  x7  +  At  *  cos  x 


*  •  *  I  ^  t  #  *  f  '  4^4  •  k  *  ^  ^  ■  •  M  *  J  *  -  • 


/oVAy 


-  &2^2X2+  'n3  “  ^2^3  ^3  C°S  *4^  *3+  ^2  ^1*2*3”  Hf~ 


i-  cos  x^  sin  x 


x^  tan  x^  sin  x^  -  x^  tan  x^  cos 
x^  cos  x^  -  x^  sin  x^. 


To  determine  stability  of  the  equilibrium  position  x^*x  "X  “0» 
x^»x°-2.50,  we  linearize  the  system  near  this  point. 


The  resulting  linear  system  is 

g 

V<2  +  *3*3  +  4-  Sl"  *7  h  Ag 

V  -  *2  *2  h  +  ^  '  *2*3  +  *3  c“  x4>{3  •  CM  x7°  «6 

V  '  t,n  X7  ?3 


Our  aim  is  to  determine  the  eigenvalues  of  this  Bystem.  After  a 
tedious  calculation,  we  can  rewrite  the  system  as  follows: 

f  -  -  174.9  r  +  .126  r  +  .276  c 

'2  ^2  43  **7 

L  -  -  3.38  r  -  5438  r  +  201  r, 

3  s2  ’2  ^6 

t  =  ~  .0436  r 

Jb  ^3 

h  “  C2 


The  eigenvalues  of  the  matrix  of  this  system  are  computed  to  be: 

X,  -  -5438,  X  -  -174.9,  X,  -  -  1.6,  X,  -  0.0016; 

12  3  4 

thus  the  system  has  a  very  weak  instability,  which  in  practice  should 


amount  to  a  neutral  equilibrium. 
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Appendix  I 


Here  we  prove  the  theorem  of  Poincare  and  Dulac  which  is  restated 
here  for  the  reader's  convenience. 

Theorem.  The  system 

x  -  Ax  +  v(x)  ,  A  -  diag(  . . .  )  (Al) 

1  n 

can  be  reduced  by  a  formal  change  of  variables  x-y+...  to  the  form 
y  -  Ay  +  w(y) 

where  all  the  vector  monomials  ym£  comprising  the  series  u(y)  are 

s 

resonant  (for  definitions  see  section  3.2). 

Proof.  Performing  the  change  of  variables  x*y  +  h(y)  in  (Al) 
(here  lowest  degree  in  h(y)  is  _>  2),  we  obtain: 

y  ■  Ay  +  Ah(y)  -  h'(y)  Ay  +  v(y)  +  h.o.t.  (A2) 

Let  k>2  be  the  lowest  degree  of  the  vector  monomials  ym^s 
comprising  the  series  v(y);  in  order  to  eliminate  these  terms  we  must 
solve  for  h  the  equation 

h'(y)  Ay-Ah(y)  =  v  (y),  (A3) 

k 

where  v  (y)  retains  only  the  monomials  from  v(y)  of  degree  k. 


Our  aim  is,  therefore,  to  study  the  properties  of  the  operator 
(Lh)(y)  -  h'(y)Ay-Ah(y) 

mapping  the  set  of  polynomial  vectorfields  h(y)  of  degree  k  into 
itself.  The  complete  information  about  L  is  contained  in  the 
following 

Lemma .  The  eigenvalues  of  L  are  the  vectorfields  of  the  form 

ym£  ;  the  corresponding  eigenvalues  are  (p,m)-y  . 
s  s 

Postponing  for  a  moment  the  proof  of  the  lemma,  we  finish  the 
proof  of  the  theorem.  Using  the  lemma,  we  can  eliminate  all 
nonresonant  monomials  of  degree  k  in  the  right  hand  side  of  (A3),  by 
the  proper  choice  of  the  terms  of  degree  k  in  h(y).  Similarly,  we  can 
kill  all  nonresonant  terms  of  degree  (k+1)  without  affecting  any  terms 
of  lower  degree  -  Indeed,  the  transformation  x**y+h  (y)  does  not 

- -  k+1 

affect  terms  of  degree  Oc,  as  follows  from  (A2). 

Inductive  application  of  this  argument  proves  the  theorem. 

Proof  of  the  lemma. 

L(ym£  )  -  (ym£  )l  Ay  -  A(ym  i  )  - 


